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Abstract: We consider the entanglement entropy for the 2D Ising model at the conformal 
fixed point in the presence of interfaces. More precisely, we investigate the situation where 
the two subsystems are separated by a defect line that preserves conformal invariance. 
Using the replica trick, we compute the entanglement entropy between the two subsystems. 
We observe that the entropy, just like in the case without defects, shows a logarithmic 
scaling behavior with respect to the size of the system. Here, the prefactor of the logarithm 
depends on the strength of the defect encoded in the transmission coefficient. We also 
comment on the supersymmetric case. 
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1 Introduction and summary 

Entanglement entropy in quantum systems has been investigated in recent years in many 
different fields, ranging from quantum information theory to black hole physics. It encodes 
the information of the entanglement of a subsystem A with the rest of the system. In 
1-1-1 dimensional systems at the critical point, the vacuum entanglement entropy for a 
subsystem A which is obtained by geometrically singling out an interval is given by [1-3] 

Sa = -TTpAlogpA = ^logL, (1.1) 


- 1 - 




where pA is the reduced density matrix and L the length of the interval specifying A. 

In this paper, we consider situations where the separation of the full system is not 
merely geometric. Rather, we investigate the case where there is a physical separation given 
by a conformal interface (to which we also refer as “defect”). This is a one-dimensional 
domain wall, localized in the time direction, separating the two-dimensional space-time into 
two parts. An interface CFT naturally consists of two sub-systems - the two theories that 
are joined by the defect line. The domain wall can be fully or partially transmissive, such 
that the quantum field theories living on the two sides are related nontrivially across the 
defect line. A useful and interesting probe to the whole system is the quantum correlation, 
i.e. entanglement, between the two sub-systems. The subsystem A in (1.1) thus consists 
of the CFT on one side of the interface. 

In the present paper, our discussion focusses on the two-dimensional Ising model, 
where defects have been analysed by integrability [4, 5] and conformal field theory [6, 7] 
techniques. There are altogether three classes of defects preserving conformal invariance. 
Two of them have a simple description in terms of a square lattice model. At the position 
of the interface, the couplings between the spins are different than in the bulk of the lattice. 
In formulas, the energy-to-temperature ratio is given by 

^ = - (AifTijfJi+ij -h (1 - b)Ki ^ (1.2) 

3 

where Ujj = ±1 are the spin variables, and sinh(2Ari) sinh( 2 Ar 2 ) = 1 so that the bulk 
theory is critical. Along the (vertical) interface, couplings are rescaled by the factor b that 
parametrizes deformations of the interface. In the special case b = 1 the situation reduces 
to the case without any defect, on the other hand, for b = Too or 6 = 0, one obtains two 
isolated subsystems, separated by a totally reflective defect. One furthermore distinguishes 
between ferromagnetic interfaces for which the parameter b takes values b G (0, oo) and 
anti-ferromagnetic interfaces which are parametrized by 6 G (—oo,0). 

In a spin chain interpretation, the defect sits on a particular link of the spin chain, and 
we consider the entanglement entropy of the subsystems located left and right of the defect 
link. When the system propagates in time, the defect link sweeps out a one-dimensional 
line in two-dimensional space-time, which is the defect line of the conformal held theory. 

One now expects the entanglement entropy between two subsystems to reduce to (1.1) 
in the totally transmissive case at 6 = 1 and to vanish in the totally rehective case. For 
generic b, the entanglement entropy will depend on b, as this parameter determines the 
“strength” of the defect. 

Indeed, we will show that the entanglement entropy is given by 

S = a{T)logL + C, (1.3) 

where C is a constant independent of L and T parametrizes the transmission of the defect. 
For the defects (1.2) 7~ can be expressed in terms of b and the constant C vanishes. The 
formula for the case of the remaining class of defects, which is not described by (1.2), is 
similar. This class can be obtained by appyling order-disorder duality to one side of the 
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interface, which does not change the transmissivity and also not the prefactor On 

the other hand, it does shift the constant C. 

Entropy formulas of the form (1.3) have appeared in several different contexts in 2D 
GET. Eirst of all, the entropy through a conformal interface was considered in the example 
of the free boson compactified on a circle in [8] with a result precisely of the form (1.3). 
Indeed, our analysis in sections 3, 4 ,5 is similar to their computations. 

The entanglement entropy for subsystems separated by a defect in the Ising model as 
well as other fermionic chains was studied before with different methods. Numerical results 
were presented in [9], subsequently an analytical analysis appeared in [10]. In particular, 
the form of the entanglement entropy (1.3) was derived using the spectrum of the reduced 
density matrix in the lattice model. The paper [10] initiated a series of following papers 
addressing related topics, see e.g. [11-15]. 

In the present paper, we investigate the defects of the Ising model from the point of 
view of conformal field theory. Here, one associates one class of defect lines to each primary 
of the Ising model. Within each class, the elements differ by marginal perturbations and 
as a result also by their transmission and we verify (1.3) by a conformal field theory 
computation. The different classes differ physically by properties such as their ^-factor and 
RR-charge. From the point of view of the entanglement entropy, this changes the constant 
contribution in (1.3). The constant shifts are particularly interesting in the supersymmetric 
case, which we analyze by combining the Ising-model results with those of [8] . 

Constant shifts in the entanglement entropy were also observed in [16-18]. In those 
papers, the shift in the entanglement entropy of excited states, rather than the vacuum, 
was determined. The excited states were obtained by acting with local operators on a 
GET vacuum. The physical difference to our situation is that the defects we consider 
extend in one dimension, hence are not local operators. In the case of the papers [16- 
18] the logarithmic term remains the same (compared to the situation of the vacuum), 
whereas the constant term gets shifted. In the case of rational conformal field theories, 
the constant shifts have an interpretation in terms of quantum dimensions. In [19, 20] the 
entanglement entropy was considered for conformal field theories with boundary. These 
systems are related to ours by the folding trick, where one folds along the defect line to 
obtain a tensor product theory with a boundary, see figure 1. However, the division in 
subsystems is different in their case, as they consider the division of the system into left 
and rightmovers and compute the left-right entropy. 

This paper is organized as follows: In section 2 we review the construction of conformal 
interfaces for the Ising model. Using the folding trick, conformal interfaces can be mapped 
to boundary conditions for a free boson on a circle orbifold. Hence, they are given in 
terms of DO branes (ferromagnetic and anti-ferromagnetic) and D1 branes (order-disorder) 
parametrized by their position and Wilson lines, respectively. While this description offers 
an intuitive interpretation of the possible interfaces, it is less useful for calculations in our 
context. We hence go back to a formulation in terms of a GSO projected free fermion 
theory. The ferromagnetic and anti-ferromagnetic interfaces are then given by interfaces 
charged under RR-charge whereas the order-disorder interface is a neutral interface. 

In section 3 we explain the basics of how to compute the entanglement entropy as a 
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derivative of a partition function involving defects. Sections 4 and 5 contain the concrete 
calculations for the Ising model, in particular (j{T) and C in equation (1.3) is derived for 
all classes of interfaces of the Ising model. In section 6 we comment on the supersymmetric 
model, combining our results with those of [8]. We show that (t{T) simplifies due to can¬ 
cellations between oscillator modes. We also express the result in a way that is suggestive 
for generalizations to higher dimensional tori and possibly other supersymmetric models. 
Finally, in section 7 we draw some conclusions and point out open problems. 

2 Conformal interfaces of the Ising model 
2.1 Interfaces and boundary conditions 

A convenient description of defects in the Ising model arises, when we employ the folding 
trick. Here, as illustrated in figure 1, an Ising model interface is mapped to a boundary 
condition of the tensor product of two Ising models. The latter is well known to be 
equivalent to a Z 2 orbifold of a free boson compactified on a circle of radius r = 1 [6, 7]. 


<—)• 


interface CFT boundary CFT 


Ising 

Ising 

CFTl 

CFT2 



Figure 1. The folding trick transforms interfaces (red) of the critical Ising model to boundary 
conditions of the c = 1 Z 2 -orbifold theory. 


The boundary conditions of the orbifold theory come in two continuous families, [6, 7, 

21 ]: 

• Dirichlet conditions with G [0,7r] , 

• Neumann conditions \N,4>'f) with (p S [0)'^/2] • 


In string theory language, 4> is the position of a DO-brane on the circle with a Z 2 
identification, whereas p is the Wilson line on a Dl-brane which belongs to the position of 
the dual DO-brane on the dual circle (of radius r = 1/2). Unfolding converts the boundary 
states \B)) of (Ising)^ to interfaces of the Ising model. 

For Dirichlet interfaces one can relate (p and the parameter b of the interface model 
(1.2) as in [6, 7] by comparing the CFT spectrum with the exact diagonalization of the 
transfer matrix [22]: 


tan(i/) — 7r/4) = 


sinh(iiri(l — b )) 


i ->■ COt { p ) = 


tanh(677i) 


( 2 . 1 ) 


sinh(Ari(l + h)) tanh(i^i) 

A special case is (/ = 7r/4 corresponding to 6 = 1 which means there is no interface. 
Hence the interface operator is given by the identity operator. Another special case is 
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(/> = 37r/4 which belongs to 6 = — 1. This operator belongs to the Z 2 -symmetry of the Ising 
model. 

At the special values (/> = 0,7r/2 and tt, corresponding to 6 = oo, 0 and —oo, respectively, 
the interfaces reduce to separate boundary conditions for the two Ising models given by 

(++)©(-), iff) and (H—)©(—h), (2.2) 

where / dennote the three conformal boundary conditions of the Ising model, namely 
spin-up, spin-down and free [23]. 

For Neumann interfaces, or order-disorder interfaces, the relation between (p and b is 
similar but with (p G [0,7r/2] and thus 6 > 0 (see e.g. in [24]). Again, for the special value 
(p = 7r/4 we have 6=1. This means that this Neumann interface is topological. On the 
other hand, at the values <^ = 0,7r/2 the interfaces reduces to separate boundary conditions 

(+/) © (-/) and (/+) © (/-). (2.3) 

Other two interesting quantities that characterise all conformal interfaces are the reflec¬ 
tion coefficient TZ and the transmission coefficient T which are given by 2-point functions 
of the energy momentum tensor as follows [21] 

_ (TiTi + T 2 f 2 ) ^ (TiTa + ^ 2 ^ 1 ) 

((Ti + f2){fi + T2)) ’ ((Ti + r 2 )(ri + T2)) ■ ^ 

where Ti, Tf are the components of the energy momentum tensor at the point z and T 2 , T 2 
are evaluated at the corresponding point reflected at the interface. For the interfaces we 
considered previously the reflection and transmission coefficients are given by 

cos^(2(/)) Dirichlet 
cos^(2(/>) Neumann ’ 

It is easy to see that TZ + T =1. Note that for topological Dirichlet interfaces, where 
(p = 7r/4 or there is no reflection, namely 7?. = 0. On the other hand, for cp = n7r/2 

the reflection coefficient is 7^ = 1, and thus the interface reduces to a totally reflecting 
boundary conditions. For Neumann interfaces the statements are alike. 

2.2 Free fermion description 

While the description in terms of the free boson provides an overview over the possible 
interfaces, to construct the explicit interface operator one needs to undo the folding. This 
is best done in the language of free fermions. Recall that the Ising model can be regarded 
as a system of a free real Majorana fermion, where modular invariance is achieved by 
a projection on even fermion number (where the fermion number is the sum of left and 
right fermion number). In a free fermion theory one distinguishes between the NS-sector 
and the R-sector. In the NS sector, the fermions 'ip,'ip (denoting left and rightmovers) 
are half integer moded and there is a non-degenerate ground state. In the R-sector the 
fermions are integer moded and the ground state degenerates. The Ising model has three 



and T = 


sin^(2(/i) Dirichlet 
sin^(20) Neumann 


(2.5) 
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primary fields with respect to the Virasoro algebra, 1 ,ct, e of left-right conformal dimen¬ 
sions (0, 0), (1/16,1/16), (1/2,1/2). In terms of the free fermion (0,0) is the NS-vacuum, 
(1/2,1/2) the first excited state of the NS-sector and (1/16,1/16) a R-ground state (after 
the degeneracy of the ground states has been lifted by the GSO-projection). 

Having an interface between two 2D free fermion conformal field theories as on the left 
of figure 1, its interface operator has the general form [24, 25] 


/i.2(0) = n Il2{0)ll2{0) = , (2.6) 

n>0 

where we have split the operator into two factors; if 2 ( 0 ) is a map of the ground states 
of the free fermion theory, whereas 1^2 contains the higher oscillator modes. The latter 
can be factorized further; in 1^2 only the ±nth modes of the fermion field appear pairwise, 
such that all If 2 commute. It is given by 


Il,2 = exp (-iV'inC’llV'-n + 'Ip-n0l2'lfn + 211^1 + ifjl022'4’l) , (2.7) 

where the ^'^e the modes of CFT1/CFT2 which are acting from the left/right on I^ 
~ the ground state operators. The matrix O G 0(2) specifies the interface and can be 
given in terms of a boost matrix A G 0(1,1) which guarantees that the interface preserves 
conformal invariance, see [24, 25] for more details. Their exact relation is given by 


0(A) 


( ^12^22 

\ ^22 


An — AioAoo^A 


12 ^ 4-22 ^'■21 
-^22^ A21 


( 2 . 8 ) 


The matrices with detA = -|-1 correspond to Dirichlet boundary conditions in the 
orbifold theory whereas detA = — 1 corresponds to the Neumann boundary conditions. For 
detA = -|-1 the relation gives 


A 

and for detA = — 1 


( coshy sinhy 
sinhy coshy 


GG O 


( cos(2(/>) sin(2iji) \ 
ysin(2(/)) — cos(2 ^)j ’ 


(2.9) 


^ /coshy — sinhy\ ^ q I cos{2(p) sin(2(/>) 

Isinhy — coshyy \—sm{2(f>) cos(2(^) 

Indeed, cj) and (j) precisely correspond to the parameters describing the DO and D1 brane 
moduli space. From now on we omit the tildes. Then we can write cos(2y)) = tanhy yy 
= cot (j) in both cases. To obtain the interface operators of the Ising model from those of 
the free fermion theory, one still has to GSO-project on total even fermion number. This 
requires taking linear combinations of the free fermion interfaces. In the Ising model the 
type-O-GSO projection allows us to distinguish three cases: The interface operators for 
det A = 1 that carry either positive or negative RR charge and can be written as ^ . 

^ See [24] for more details on the construction 
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(2.11) 


/±(A) = 2 ± I^{A)) + (A ^ -A), 

and these for det A = — 1 which are the neutral operators 

/°-(A) = ^J^S(A) + (A^-A). (2.12) 

v2 

The operators and act on the Neveu-Schwarz and Ramond sector of the free 
fermion theory, respectively. They are given by 

/NS(A)= Yl r(A)/0’NS^ with /°’''S = |0 )nsns(0|, 
neN-| 

and 

/R(A) = r{A)I°’^, 

nGN 

with 1°’’^ = \/| sin(2(/>)| (^|+)rr(+| + |-)rr(-|)5(A) 

= \/2(^cos((/>) |+)rr(+| +sin((?i) |-)rr(-|^ . 

Here, |±)r denote R-ground states and S is the spinor representation of 0(1,1). 

3 How to derive the entanglement entropy 

In this section we briefly want to show the derivation of formulas we use in the following sec¬ 
tions. It is an adaptation of the procedure used to derive the entanglement entropy through 
interfaces in the free boson theory. Thus, for a more detailed derivation we recommend [8], 
but also [2]. 

Our goal is to derive the (ground state) entanglement entropy between two 2D GSO- 
projected free fermionic CFTs connected by the interfaces (2.11) or (2.12). We formally 
define CFTl to live on a half complex plane Retc > 0 and CFT2 on Rerc < 0, respectively. 
The interface then lies on the imaginary axis Hew = 0. The EE is defined by the von 
Neumann entropy of the reduced density matrix for the ground state pi = Tr2|0)(0| as (see 
e.g. [1, 26]) 


(2.13) 


(2.14) 


5 =-Tripilogpi =-dicTripf . (3.1) 

The trace of the AT-th power of the reduced density matrix is also given by the partition 
function on a iF-sheeted Riemann surface TZk with a branch cut along the positive real 
axis [2, 8] 
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(3.2) 


Tripf 


Z{K) 

Z(l)^ 


The iC-sheet construction is illustrated on the left of figure 2. This procedure is a 
version of the so-called replica trick. By the use of (3.2) the entanglement entropy can be 
written as 


S = {l-dK)\ogZ{K)\K^l. 


(3.3) 


Im w 

Interface 
branch cut 

Re w 


Im z 



Figure 2. Sketch of the if-sheet Riemann surface TZk that one needs for the replica trick. After 
imposing a IR citoff e, an UV cutoff L, and suitable variable transformation, TZk looks like on the 
right. 

To evaluate the partition function Z{K) we change coordinates lo z = logtc and 
introduce the cutoffs \w\ = e and |rc| = L. The ii'-sheet TZk then looks as on the left part 
in figure 2. As in [8] we impose periodic boundary conditions in Re z and choose e = j; for 
simplicity. We then end up with the torus partition function with 2K interfaces inserted 
which is given by 


Z(K) = TVi (li,2 l2,i • • • /2,i 

= Tn (/i ,2 42 (3.4) 

with S = 7r^/(logL) and / 2 ,i = 42- 

4 Derivation of the partition function 

In the following we explicitly derive the partition function (3.4) for the interface operators 
introduced in section 2.2. We start with a single NS operator which is the simplest case. 
Step by step we show how to derive Z{K) for the more complicated R, neutral, and charged 
operator. 
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4.1 The partition function for a single NS operator I = /^^(A) 

In the NS sector of the free fermion theory we can formally write the Hilbert space % of 
the theory as a tensor product ®nHn where Hn = span 110), V'-n^-n|0), V'-n|0), ^-n|0)}. 
Then each I'^{0) as in (2.7) has a matrix representation in T-Ln given by 

/ 1 -iC>22 0 0 \ 

_ -i Oil - det O 0 0 

0 0 c>i2 0 ’ ^ ^ 

y 0 0 0 O 21 J 

and we can write In this notation the propagator is given by = P = 

®nP^ with 


pn 


diag(l,e 


— 26n ^—Sn ^—5n\ 


(4.2) 


Using the above notation with CFTl = CFT2 the partition function (3.4) of the K- 
sheet can be written as 


z{K)= n 


(4.3) 

nSN-i 



= n 

(A^l + A^2 + A^3 + Xn,4:) ’ 

(4.4) 


nSN-i 


where Xn,i, i = 1, 4, are the eigenvalues of D"‘ = [/"■ P” (/"')tP’^]. 


4.1.1 Explicit calculation 

We have two distinguishable interfaces: the Dirichlet-interface with detA = 1 and the 
Neumann-interface with detA = —1. However, in both cases the matrices Zl"' are similar 
and their eigenvalues are given by 


An,i = e = e ^cosh(2n(5) + cos^(2</)) + cosh(re<5)-\/2 cosh(2n(5) + 2cos(4(/))^ , 

An ,2 = ^cosh(2n(5) + cos^(2</)) — cosh(re<5)-\/2 cosh(2n(5) + 2cos(4(/))^ , 

An ,3 = sin^(2())) = An,4 , 

so that the partition function of the AT-sheet is given by 

Z{K)= n + + . (4.5) 

nSN-i 

At this stage one could proceed further by directly using formula (3.3) on the latter 
result for Z{K). Through the logarithm the infinite product simplifies to a sum. Taking 
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the derivative w.r.t. K in every summand it is then easy to write down a result for the 
entanglement entropy by means of an infinite sum. One could then evaluate the sum - 
and thus the entanglement entropy - numerically for every <5 > 0 and 4> up to arbitrary 
accuracy. However, we are mainly interested in small 5 - which means large L ~ behaviour 
of the entanglement entropy, since L is introduced as a UV cutoff. In this limit we can 
derive the EE analytically by proceeding as in the following. 

Eor odd K the partition function (4.5) can be written as 

Z{K)= (2cos2(i/fc) - 1+ cosh(2n5))V (4.6) 

neN-| \k=l / 

with Uk = arcsin (sin(2(/>)| sin(^)l). Eor even K we have to add sin^^(2(/>) to 

every factor in (4.6). Additionally, we state that the fraction 0[0,0](r, 2 )/r/(r) of the well 
known 0-function and r/-functions as defined in (A.l) and (A.3) can be written as 

= elf n (2 cos^ (vrz) - 1 + cos (2n7rr)) . (4.7) 

^ ^ neN-| 

Thus we can conclude that the AT-sheet partition function for odd K can be expressed 
as 


, . KS 

Z{K) = e“i2- 


n^Mini 
M »(f) 


Using the behaviour of r] and 9 under S-transformations we can write 


(4.8) 


Z{K) = 


KS 
■ 12 


K 

H' 

k=l 






Z{K) = ew e- 


s 


1 + e' 


(4.9) 


where (p{K) = ^ki constant in 5. 

Eor even K, Z{K) can not be given in terms of 0 and r/ as above. One might wonder 
if we really can use (4.9) to calculate the entanglement entropy although it is just valid for 
odd K. In Appendix B we actually show why it really suffices to consider (4.9). 


4.2 The partition function for a single R-operator / = /^(A) 

In the Ramond sector the modes are integer and the zero-mode map is slightly more 
difficult, = \/2f cos(</>) |+)rr(+| -|-sin((/)) |— )rr(— l). The latter allows us to write 


= cos(^) + sin((/>) , with /^ = \/2 |±)(±| I” , (4-10) 

neN 

where /5'(A) • /^(A') vanishes. Proceeding similar to the case of NS-operators one gets 
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(4.11) 


Z{K) = 2^ {cos{(f)f^ + (A^i + A^2 + -^^3 + > 

nSN 

where again An,i, i = 1, ••, 4, are the eigenvalues of [/" P” (/”)'^P”]. 

4.2.1 Explicit calculation 

The eigenvalues for the R-interface are similar to the eigenvalues for the NS interface but 
with n G N. Thus, for odd K we can write 

Z{K) = 2^ (cos((/))^'^ + sin((/))^^) f (2cos^(t'fc) — 1 + cosh(2(5n))^ , 

fc=l VnSN / 

(4.12) 

where again = arcsin (sin(2(/>)| sin(^)|). This time the latter is given in terms of 
0](r, z)/? 7 (r) because of n being integer. One important difference between the 9- 
function we use here and the 0-function used in the case of the NS-interface is that there 
appears an additional factor of 2cos(7r2:). One can show that 

K 

cos{uk) = cos{(j))‘^^ + 

k=l 

for odd K so that in this case the partition function reduces to 


K 


Z{K) = e 6 (cos(0)^^-|-sin((/))^'^) 


TT«l5.0l(f.?) 


.1 ti " (¥) 


= e 6 


(f,f) 


(4.13) 


fc=i ^(t) 

The same steps as for the NS interface now lead us to 


(4.14) 


in the limit 5 <C 1. 

In the Ramond sector only the Dirichlet interfaces have non-trivial components, so 
we do not consider Neumann boundary conditions, although they would not make any 
difference for Z{K). 


4.3 The partition function for the neutral interface operator 

The neutral interface operator is given by /’^■(A) = ^ (l^®(A) -|- I^®(—A)) with Neumann 
boundary conditions, i.e. det A = —1. Some simple algebra leads to 
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+ iNS(_a) e-SH jNS(a) e-5H ^ I^^{-A) e-^^ 


= 2 (iNS(a) /NS(a) ( 4 .I 5 ) 

+ 2(/NS(A)e-'5^/NS(_A)e-5/^) 

= 2 (L>+ + L)_) , 

with D± given by the tensor product of the matrices 


Dl 


( 1 + n 

—i{e~‘^^"‘ + e 
0 

V 0 


i{e 


—4Sn _|_ ^—6Sn 


)Vn 


-45n)^ g-45n_^g-25n^ 

0 
0 


0 
0 

Ae-2<5n J- 
0 


±e 


0 

0 

0 

—25n 


(4.16) 


'^7 


where we here used the reflection coefficient TZ = cos(2(/))^ and the transmission coefficient 
T = sin(2(/))^ as introduced in (2.4). We can see that both and can be diagonalized 
simultaneously. With straight forward linear algebra we can now calculate the partition 
function for the X-sheet with a neutral interface insertion. It is given by 


Z{K) = Tr(A>+ + Z1_)^ 

= 2^-1 n e-'''"'(+2sin2^(2<()) + (p+)^ + (p-)^) + 

neN-| (4-17) 

+ 2^-1 H e-^^-^{-2sm^^{2cP) + {pif + {p-f). 

neN-| 

The important factor of 2^~^ can be understood with the following simpler example: 
Consider the matrices M±{x) = (8>idiag(x, ±x). It is now easy to convince oneself that 

[M+(x) + M_(x)] • [M+( 2 /) + M_(y)] = 2[M+(xy) + M_(xy)], 

which allows us to directly conclude that [M+(x) + M_(x)]'^ = 2^“^[M+(x^) + M_(x'^)]. 
The generalization to Z1+ and D- is straight forward. 

The first summand in (4.17) is the same as for the single NS-interface. In a very similar 
way as in section 4, the second summand can also be written in terms of 0[0,0]/r/. We here 
only state the result in the limit <5 <C 1: 


nK—^ / _‘pCP 1 x{k) \ 

Z{K) = 2^ 125 (e ^ +e s j 


(4.18) 


where xi^^) = J2k=i = arcsin(sin(2(/))| cos(^)|) 
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4.4 The partition function for the charged interface operator 

The charged interface operator is given by /^(A) = ^ (/^®(A) ± /^(A)) + (A —)• —A) with 
Dirichlet boundary conditions, i.e. det A = 1. With similar considerations as for the neutral 
interface we can write 


27^(A) /^(A)t = {D+ + £>_)© 2 cos^{(^) (^77^+ + 7?^+) © 

©2 sin2((/,) [dI- +D^-'^ , 


(4.19) 


where corresponds to the vacuum |+) and D^~ corresponds to |—) in a similar way 
to the single Ramond interface. There is no difference between the positively and the 
negatively charged interface. In matrix representation all the D’s are tensor products of 
matrices similar to (4.16) but with integers for the Ramond operators. Thus the partition 
function of the K-sheet with a charged interface can be written as 


Z{K) 


^Tr 

2K 


{D+ + D- f ® 2^ cos{cl)f^ ( 


7jR+ ^ ^R+\ 


K 


© 2^ sin((/)) 


2K 


T^R jjn 


R- 


K 1 


= l( U e-™(+2sin2^(2<()) + (p+)^ + (p-)^) + 

'^neN-i 

+ n {-2sm‘^^{2(l)) + {p+)^+ {p-f) + 

neN-| 


(4.20) 


+ 2^ (cos(()))2^ + sin(0)2^) J] (+2 sin2^(2</.) + {p+f + (p-f) + 

nSN 

+ 2^ {cos{(l)f^ + sm{cj)f^) (-2sin2-^(2<^) + (p+)^ + (p")^) 

nSN 

Using the same logic as in the previous sections, the partition function for odd K 
reduces to 


„ ,, , X(.^) 

Z{K) = e w e-^ + f{K) 


(4.21) 


in the limit <5 <C 1. The functions 4> and y are given as before and f{K) is given by 

. / K 


1 


/(A') = -\^l + + sm{0)“-) n 


/I 


(4.22) 


5 Derivation of the entanglement entropy 

Before we explicitly derive the entanglement entropy we want to show that - for <5 <C 1 - 
it is the same in all previous cases up to an additional log 2 term for the neutral interface. 
Our formula of choice (3.3) is 
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S = {l-dK)logZ{K)\K^^, (5.1) 

for which it is easy to check that any overall factor in Z{K) with C constant in K 
does not contribnte to the entanglement entropy. 

At first we want to write down the preliminary result for the single NS-interface where 
Z{K) is given by (4.9) so that the entanglement entropy can be written as 

5NS = (_^(1) +5x^(1))!. (5.2) 

Next we consider the single R-interface where the partition function is given by (4.14). 
Its EE is simply the same as for the NS-interface 




(5.3) 


Now we want to derive the EE for the neutral interface. Inserting its partition function 
(4.18) in (5.1) gives 


S"- = log + 

V J 6 e-v?(i)/(5 _|_ g-x(i)/5 

= (-(/9(1) dK(p{l)) ^ - log 2 
= 5^^ -log2L 


log 2 


(5.4) 


where we can simplify to the second line because we are in the limit d <C 1 and because 

< 4 >( 1 ) = 0 > x(l) = 

At last we want to derive the EE for the charged interface operator where Z{K) is 
given by (4.21). As in the case of the neutral interface operator every term with a factor 
g-x(i)/5 can be neglected. Consequently also f{K) in (4.21) has no contribution to the 
entanglement entropy when we are in the limit <5 <C 1. It again simply reduces to the EE 
for the single NS-interface: 


S± ^ ^NS 


(5.5) 


5.1 Explicit derivation of 

To derive the entanglement entropy explicitly we have to calculate Therefore we 

proceed similar as in [8] and write ip{K) = Ylk=i / (f)) which can be written 

as a Taylor series around k/K = 0 and further massaged as 


K oo 


if{K) = EE/» 


k=l m=0 

oo 


K 


OO « K 

/ ^ ]^m / j 

m=0 k=l 


fm -|- 1) — Bm+1 

/ ^ iX™ 


m=0 


m -|- 1 


(5.6) 


- 14 - 







where Bn{x),Bn are the Bernoulli polynomials and Bernoulli numbers, respectively, as 
given in Appendix A.3. Its derivative in the limit A —)• 1 is then given by 


dK'^{K)\K^l = , dKBm+l{K + l) \k^I — 1 (-Bm+l(2) — Bm+l) 

Tn ^ ^ 

=dK{Bm+i{K) + {m+l)K^) =(m+l) 

= 1 dKBm+l{K)\K^l + fmm- frnm 

m+1 

m 

= y^^^dKB^+i{K)\K^i. 
m+1 

m 

(5.7) 


At this stage we use the formula (A.8) to obtain 


1 r 

dK^{K)\K^i = /(O) + 2 /( 0 ) + J 


if (it) - if {-it) 


o2TTt 


dt. 


10 1 - e' 

Both /(O) and f'{0) vanish, so that the entanglement entropy is given by 


(5.8) 


SNS ^ 1 ifiif-ifj-if ^^ 


S JO 
Tr‘^a{s) 


27Tt 


1 - e 


= cris)log{L) 


(5.9) 


where we defined s = | sin(2(/>)| = ff and 


o roo arcsinh(s sinh(7rt)) (coth(7rt) — l) s cosh(7rt) 

a{s) = - / - ^ , A - 

^•^0 + s2 sinh2(7rt) 

=^r“(\/^Aij+ 

For the last step we substituted u = arcsinh(s sinh(7rt)). 

5.2 Comments on the result and special cases 

The result (5.9) shows that the EE through a single NS interface has a logarithmic scaling 
with respect to the (large) size L of the system. Up to an additional term “log 2” for 
the neutral interface operator this is exactly the behaviour of the EE through a conformal 
interface in the 2D Ising model, too. The interface itself affects the EE mainly through 
the factor cr(s) which is given in integral form in (5.10). The square of the “variable” 
s = I sin(2(/))| just is the transmission coefficient of the interface T, but can be given in 
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terms of the scaling factor b and the coupling constant Ki of the Ising model as in (2.1), 
too. In figure 3, we show the explicit dependence of the factor a on the transmission 
coefficient T. 



T 

Figure 3. The prefactor cr as a function of the transmission coefficient T. As expected, it vanishes 
for totally reflecting interfaces, where T = 0, and reproduces the familiar result a{T = 1) = 1/6 = 
c/3 for topological interfaces in the Ising model. 


There are two special cases to mention. First, we want to consider topological interfaces 
where the transmission coefficient is maximal, T = 1. This case includes the identity in the 
free fermionic CFT and in the Ising model which is why it should reproduce the universal 
scaling of the entanglement entropy without interfaces. In fact, with (t( 1) = 1/6 = c/3 we 
obtain the right result from [1]. However, not just the universal scaling is correct, but also 
the sub-leading term, i.e. contributions constant in L, are right. This can be checked by 
directly evaluating the torus partition function in the theory without interface, which is 
e.g. for the Ising model given by 


1 0[|,|](r,O) 0[i,O](r,O)\ 


(5.11) 


with T = i^. 

TT 

Secondly, in the limit of totally reflecting interface, when T —t 0, one can show that 
a{T) = 0{T\ogVT), so that for vanishing transmittance also the entanglement entropy 
vanishes. This fits the fact that all the oscillator parts of the two CFTs are decoupling for 
an interface with T —)• 0. As shown in figure 3, cr{T) increases monotonically between these 
two extremal cases. The latter observation supports the intuitive guess that entanglement 
changes according to the transmittance of the interface. The lower the transmittance the 
lower the strength of interaction between the two CFTs connected by the interface. 

We also want to discuss the sub-leading term although it vanishes but for the neutral 
interface operator. An important contribution to that term for the Ramond and the charged 
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interface comes from the ground state map.^ In our cases the ground state map is rather 
easy, it just separates the Hilbert space of the CFT in direct sums. Let us consider a 
Hilbert space T-L = I-Lq © Hi and two density operators po and pi acting on the respective 
Hilbert space, such that 


p = apo©/3pi, (5.12) 

where a and /3 are real numbers. It is then easy to see that the von Neumann entropy on 
the full Hilbert space Sp can be given in terms of the entropies on Hq and Hi: 

Sp = aSpo + f3Sp^ - {a log a + f3 log f3). (5.13) 

As an example let us consider the single Ramond interface. Because of the ground 
state map the (reduced) density matrix can be written as p = cos^(/> + sin^c/) /?_ and since 

Sp^ = Sp_ the von Neumann entropy reads 

= Sp = Sp^ — {cos‘^(j) logcos^c)) + slv^cf) logsin^(/>). (5.14) 

However, Sp^ itself has an additional sub-leading term that exactly cancels the contri¬ 
bution of the ground state map. 

Another contribution to the sub-leading term comes from the GSO-projection which 
separates the Hilbert space of the free fermion in a direct sum H = Ho ®Hi, graded by 
the fermion number. As an example consider the interface operator /^®(A) on the full NS 
Hilbert space. Then the projected interfaces on Ho and Hi are 

/± = i(;”®(A)±;"®(-A)), (5.15) 

where /+ is the previously called neutral interface operator. A density operator is given 
by /!/, such that p = \{p+ + P-) and thus 

= Sp = ]^{S+ + S.) + log2, (5.16) 

where in the limit of large L it can be shown that 5+ = 5_, so that we get = 5“' ©log 2 
or 


-log2. (5.17) 

It is noteworthy that the sub-leading term does not depend on the most significant 
property of the interfaces, namely the transmittance 7~, in all our cases. A similar result is 
known for entanglement entropy through interfaces in the free boson theory [8] . There the 
sub-leading term only depends on the winding numbers ki , /c 2 on both sides of the interface 
and is simply given by log |/ciA; 2 |”. 

^Especially for bosonic interfaces this is often called the lattice part of the interface operator. This name 
comes from the change of the lattice structure of the torus by the ground state mapping. 
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As a final comment in this section we want to state that one can also express the 
pre-factor a in terms of the Logarithm and Dilogarithm Li 2 in a similar way as in [8]. It 
then looks like 

cr(s) = ^ ^ - ^((s + l)log(s + l)logs + (s - l)Li 2 (l - s) + (s + l)Li 2 (-s)) . (5.18) 

This indeed agrees with formulas (22-26) of [10].^ 


6 Supersymmetric interfaces 


Let us now consider a situation with = 1 supersymmetry by adding a free boson to 
the theory of a free fermion. We can combine our results with those of [8] to obtain the 
entanglement entropy through a supersymmetric interface. Compatibility of the interface 
with A" = 1 supersymmetry requires that the interface intertwines the supercurrents G^, 
of the two theories 

(G^ - i4Glr)h,2 = hMGl - irilGl,) , (6.1) 

where rj, rjg, rjl = ±1. The signs rjg specify the preserved SUSY of the two theories and the 
GSO projection requires to sum over both choices for r] in the final step. As explained in 
[25] the choices of sign can be absorbed in the gluing matrix A G 0(1,1) for the fermions, 
such that specific entries in the gluing matrix for bosons and fermions can differ by signs. 
Since the entanglement entropy does not depend on these choices, we simply assume that 
bosons and fermions are glued by the same matrix A (or equivalently O) that we used in 
the current paper, implying that the preserved SUSY is the same on the two sides of the 
interface. 

The implementation of the GSO projection has been discussed in section 2.2 and can 
be taken over to the supersymmetric situation. 

The full interface operator of the supersymmetric theory can be written as a tensor 
product of a bosonic and a fermionic piece: 

Ii,2(A, ki, k2, if) = /i,T(A, ki,k2, v?) 0 /q2“(A) • (6.2) 


Interfaces of the theory of a free boson compactified on a circle were considered in [8, 27]. 
They depend on two integers, ki and k 2 , specifying topological winding numbers, as well 
as two continuous moduli, (pi, ^ 2 - The origin of these parameters is easiest to understand 
in the “folded” picture, where the interface is mapped to a D-brane on a torus. In the 
simplest case of a one-dimensional brane wrapping the torus x S^, the integers can be 
understood as winding numbers of the Dl-brane around the two 1-cycles. The continuous 
moduli correspond in this picture to position and Wilson line. The gluing matrix A is 
restricted by the torus geometry 


tanhy 


(fcli^l)" - {k2R2? 

{klRlf + {k2R2?' 


e 


®In our result there appears an additional factor of 2 since we identify the IR and the UC cutoff via 
= 1/L. 
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We choose the fermionic gluing matrix to be the same as the bosonic one. The entanglement 
entropy has been computed in [ 8 ] with the result 

3^^°^ = ais)^^nogL -log\kik 2 \ , (6.3) 


where 




s 

2 



1 \ du. 


(6.4) 


The Hilbert space of the supersymmetric theory is the tensor product of the bosonic and 
fermionic Hilbert spaces. The partition functions Z{K) hence takes the product form 
Z{K) = Z^°^{K) ■ Z^^™{K). Due to the logarithm in (5.1) the entanglement entropy can 
then be written as 


^SUSY 


^bos _|_ ^ferm 


(6.5) 


where in the fully supersymmetric model is given in (6.3) and 5^®™ is = 5^.^ 
Explicitly, 

-SSUSY = |iog_L_log|/i;iA:2| ■ (6.6) 

The prefactor (t(s) of the logarithmic term simplifies significantly. Here, contributions 
from the oscillators of the bosonic and fermionic part of the system cancel out in the limit 
(5 —)• 0, such that only the term s/2 remains. This is similar to the computations in [25], 
where the limit of two parallel interfaces approaching each other was considered. Note 
that the constant contribution log |/i;iA: 2 | has a topological interpretation: The winding and 
momentum modes of the compactified boson are quantized and part of a lattice. The 
combination |A;iA: 2 | is the index indK of the sublattice of windings and momenta to which 
the interface couples. Here, A G 0((i, djQ) is the gluing matrix for integer charge vectors 
and, as opposed to A, does not depend on the moduli. On the other hand, the quantity 
s = I sin 2(f)\ = 'JT specifies the precise geometry of the sublattice and determines the 
transmissivity, which is the same for bosons and fermions (for equal A) and also in the 
supersymmetric system. The supersymmetric entanglement can thus be rewritten as 


5 SUSY _ logL — log indA . ( 6 . 7 ) 

o 

It is very suggestive that this form of the entanglement entropy generalizes to supersym¬ 
metric torus compactifications in higher dimensions. As was shown in [25] the index of the 
sublattice is a useful quantity to characterize topological information of an interface, in 
other words, the information that does not change under deformations of the interface or 
bulk theories. In a similar way, T naturally exists for any interface and characterizes the 
transmissivity, in other words, how far away the interface is from being topological. In the 
case of higher dimensional tori, A G 0{d, d) is a 2d x 2d matrix consisting of d x d blocks 
and the transmissivity is given in terms of the determinant of the lower-right block A 22 , 

r = ldetA22l-A 

^One can also consider the GSO projection of the supersymmetric model. It has the same structure but 
with given by S'"' or S^. In the final result, there appears an additional contribution log 2 for the 

neutral interfaces, as discussed before . 
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7 Conclusions and Outlook 


In this paper we have discussed the entanglement entropy through conformal interfaces for 
the Ising model - i.e. a free fermion theory - and for supersymmetric systems. We have 
computed the prefactor (t(T) in equation (1.3) and seen explicitly how it arises purely from 
the contribution of higher oscillator modes. These largely cancel against bosonic modes in 
the model with supersymmetry. 

It would be very interesting to generalize these findings further. It is suggestive that 
also in more general systems topological data of the defect will enter the constant shift C 
in (1.3), whereas oscillator data enters the prefactor of the logarithmic term. 

The defects we investigated in this paper can be regarded as marginal perturbations 
of the topological defects of the Ising model. The latter are labelled by the primaries 
of the theory [28], in the case at hand l,cr, e. The perturbing operator is a marginal 
defect perturbation, living only on the defect. It would be interesting to consider more 
generally the entanglement entropy for initially topological defects perturbed by marginal 
and possibly also relevant defect operators. 

Another form of perturbation appears for the interfaces in the free boson theory consid¬ 
ered in [ 8 ] and in section 6 above. They come in several classes characterized by topological 
data fci, /c 2 , which cannot be changed under perturbations. However, for /ci, k 2 fixed, there 
are again interfaces related by perturbations, but this time marginal bulk perturbations 
deforming the CFT at one side of the defect.^ It would be very interesting to generalize 
this to other systems related by RG domain walls [29-32] (where the free boson “RG do¬ 
main walls” are obtained for /ci = A :2 = 1 in the above discussion) and to compute the 
entanglement entropy for them, e.g. in perturbation theory, similar to the discussion of 
the defect entropy in [31]. 

Apart from being interesting from the point of view of the physics of impurities, this 
program might also be interesting from the point of view of the physics of defects. In the dis¬ 
cussion of defects, well-known quantites are the 5 -factor and the reflectivity/transmissivity. 
The entanglement entropy might provide another useful characteristic of an interface, where 
the transmissivity enters the prefactor of the logarithmic part whereas topological features 
enter separately, namely in the constant part. This is different for the 5 -factor, where both 
topological data and oscillator data where they calculate the entanglement entropy up to 
a factor of two that appears since we identify the UV and IR cutoffenter in one factor. 


®There are also marginal defect perturbations for the free boson, which however change neither trans¬ 
missivity nor entanglement entropy. 
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A Special functions 

In the following we use q = 


A.l The Dedekind 77 -function 

The Dedekind 77 -function is dehned as 

00 

r}{T) = (1 - g”) • (A.l) 

n=l 

It behaves under T- and 5-transformations as 

77 (t 1) = eT5 r 7 (r), 77 = V-iTr]{T). (A.2) 

A.2 The 0-functions 

A most general form of the 0-functions is given by 


e [a, /3](r, z) = Y^ ^ 2 ni{n+a)iz+^) 

2 ^ 

= 77(r) q\-h, (^1 -h _g 

n=\ 

(A.3) 


Its modular transformations are given by 


9[a,l3]{T+ l,z) = e 6[a,a + (5 -]-]{t,z) , 

T 2 (A.4) 

e[a,l3]{ —,-) = \/^e2™^+“^0[/I,-a](T,z). 
r r 

A.3 Bernoulli polynomials and numbers 

The Bernoulli polynomials Bn{x) are defined by 

! -/-Tt 

, = y^^Bnjx)— with |t|<27r. (A.5) 

— 1 ni 

n=0 

The Bernoulli numbers are given by the polynomials evaluated at x = 0, namely 
Bn = B„(0). Odd Bernoulli numbers vanish. 

The sums of TTrth powers of integers can be expressed by the use of Bernoulli polyno¬ 
mials and numbers as 


_ Bm+l{N -I- 1) — Bm+l 

^ 771-1-1 

k=l 
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Using the facts that B'^{x) = nBn+i{x) and -Bn(l) = Bn for n 7 ^ 1 and -Bi(l) = 1/2 
together with the integral representation 


B 2 n = 4n(-l)^ 


t 


2n-l 


1 — e 


27Tt 


dt 


(A.7) 


one can conclude that 


d:,;Bn+l{x)\x^l = Sn,0 + “ (“O”) 


nf 


n—1 


n + 1 


1 — e 


27Vt 


dt. 


(A. 8 ) 


B The partition function for odd K is enough 

We here want to show that it suffices to consider the formula for Z{K) as given in (4.6) 
for odd K to derive the entanglement entropy. We assume that the natural analytic 
continuation of Z{K) as given in (4.5) 

ZiK) + (p-f + 2 (sin(2/.)e-^)''" = [] F„(K), (B.l) 

n>0 n>0 


gives the right result for the entanglement entropy.® For odd K the partition function is 
equivalently given by (4.6) whereas for even K we have to add 4 (sin(20)e“”'^) to every 
factor. Thus the analytic continuation of the partition function has the form 


z(K) s n 


n>0 . 


2K ^ 

I{K)A ^sin(2(/))e“”'^^ + JJ (2cos^(r'fc) — 1 + cosh(2n5)) 


k=l 


= \{Hn{K) + Gn{K), (B.2) 

n>0 


where I{K) is an analytic function interpolating between the values for odd and even K 
with 


I{K) 


0 for odd K 
1 for even K 


(B.3) 


We do not know the explicit form of I{K). However, we can show that it suffices to 
consider Gn{K) in (B.2) and that T(l) must vanish. Let us therefore derive the EE with 
the help of (3.3) for (a) and (b): 


®This really is an assumption since the continuation is not unique. 
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S={l-dK)\ogZ{K)\^^^ 

= E (“>«(G„(i) + if„(i)) - o;;(f )- ;^ ( i) (oUi) + Kwi) 

= Y,(^ogFn{l)-^{K{l) + G'^{l))\ . 

n>0 ^ J / 

In the last step we use that Hn{l) = 0 and Fn{l) = Gn(l)- In the following we show 
that also G'„(l) = F^{1). Then it really suffices to solely consider Gn{K) in (B.2) and, 
thus, we can derive the EE with the formula for odd K. In addition we then have to require 
that H'{\) = E(l)4 (sin(2i;/))e“”'^)^ = 0 which also implies that I'{1) has to vanish. 

To do so we start with rewriting 


(B.4) 


(B.5) 


Gn{K) = ]^ 2e ^I + cosh(2n(I) — sin^ sin^(20)^ 

= exp '^gnik/K) , 


\k=l 


(B.6) 


where gn(x) = log [2e (l + cosh(2n(I) — sin^ (vrrr) sin^(2(^))] can be expanded around 
X = k/K = 0, so that 


0 „W=exp( 2 :|=f;t'”') 

\rn >0 k=l J 

E Qnm ^m+l(^ 1 ) -^m+l 

1^ m + 1 

m>0 



(B.7) 


B]\r(u) and Bjsf are the Bernoulli polynomials and numbers, respectively. Proceeding 
similar as in (5.7) and the following one can calculate 


Gni^^) = exp 


9nm + 1) B^^i 

m>0 



Qnm 

m + 1 


dKBm+l{K)j , 


(B.8) 


which can be massaged further in the limit TT —?■ 1 as 
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G'{1) = exp (5n(l)) • ( 5'n(0) + ^ ^t) 


= 2e 


—2nS 


(1 + cosh(2n(5)) ( log 


2 e 


0 

—2nS 


(1 + cosh(2n5)) 


f°° 27r(coth(7rt) — 1) sin(2(/)) sinh(27rt) 
/o 1 + cosh(2n5) + 2 sin^( 2 (^) sinh^( 7 rt) 


dt 




At this point one only needs tedious algebraic deformations to show that the latter 
equals 


F^(l) = e ^log(e + log(e )p + 21 og(e sin^( 2 (/))) sin^ 

with p^ given as in section 4.1.1. 



Acknowledgments 

We like to thank Cornelius Schmidt-Colinet and Sebastian Konopka for helpful discussions 
on the topic. This work was partially supported by the Excellence Cluster Universe. 

References 

[1] P. Calabrese and J. Cardy, Entanglement entropy and conformal field theory^ J.Phys. A42 
(2009) 504005, [arXiv;0905.4013]. 

[2] C. Holzhey, F. Larsen, and F. Wilczek, Geometric and renormalized entropy in conformal 
field theory, Nuel.Phys. B424 (1994) 443-467, [hep-th/9403108]. 

[3] G. Vidal, J. Latorre, E. Rico, and A. Kitaev, Entanglement in quantum critical phenomena, 
Phys.Rev.Lett. 90 (2003) 227902, [quant-ph/0211074]. 

[4] M. Henkel, A. Patkos, and M. Schlottmann, The Ising quantum chain with defects (I). The 
exact solution, Nuclear Physics B 314 (Mar., 1989) 609-624. 

[5] D. B. Abraham, L. F. Ko, and N. M. Svrakic, Transfer matrix spectrum for the finite-width 
Ising model with adjustable boundary conditions: Exact solution. Journal of Statistical 
Physics 56 (Sept., 1989) 563-587. 

[6] M. Oshikawa and 1. Affleck, Defect lines in the Ising model and boundary states on orbifolds, 
Phys.Rev.Lett. 77 (1996) 2604-2607, [liep-th/9606177]. 

[7] M. Oshikawa and 1. Affleck, Boundary conformal field theory approach to the critical 
two-dimensional Ising model with a defect line, Nuel.Phys. B495 (1997) 533-582, 
[cond-mat/9612187]. 

[8] K. Sakai and Y. Satoh, Entanglement through conformal interfaces, JHEP 0812 (2008) 001, 
[arXiv: 0809.4548]. 


- 24 - 











[9] F. Igloi, Z. Szatmari, and Y.-C. Lin, Entanglement entropy with localized and extended 
interface defects, Phys. Rev. B 80 (July, 2009) 024405, [arXiv;0903.3740]. 

[10] V. Eisler and I. Peschel, Solution of the fermionic entanglement problem with interface 
defects, ArXiv e-prints (May, 2010) [arXiv; 1005.2144], 

[11] P. Calabrese, M. Mintchev, and E. Vicari, Entanglement entropy of quantum wire junctions. 
Journal of Physics A Mathematical General 45 (Mar., 2012) 105206, [arXiv: 1110.5713]. 

[12] I. Peschel and V. Eisler, Exact results for the entanglement across defects in critical chains. 
Journal of Physics A Mathematical General 45 (Apr., 2012) 155301, [arXiv: 1201.4104]. 

[13] V. Eisler and I. Peschel, On entanglement evolution across defects in critical chains, EPL 
(Europhysics Letters) 99 (July, 2012) 20001, [arXiv: 1205.4331]. 

[14] V. Eisler, M.-C. Chung, and I. Peschel, Entanglement in composite free-fermion systems, 
ArXiv e-prints (Mar., 2015) [arXiv: 1503.0911]. 

[15] M. Collura and P. Calabrese, Entanglement evolution across defects in critical anisotropic 
Heisenberg chains. Journal of Physics A Mathematical General 46 (May, 2013) 175001, 
[arXiv: 1302.4274]. 

[16] M. Nozaki, T. Numasawa, and T. Takayanagi, Quantum Entanglement of Local Operators in 
Conformal Eield Theories, Phys.Rev.Lett. 112 (2014) 111602, [arXiv: 1401.0539]. 

[17] M. Nozaki, Notes on Quantum Entanglement of Local Operators, JHEP 1410 (2014) 147, 
[arXiv: 1405.5875]. 

[18] S. He, T. Numasawa, T. Takayanagi, and K. Watanabe, Quantum dimension as 
entanglement entropy in two dimensional conformal field theories, Phys.Rev. D90 (2014), 
no. 4 041701, [arXiv: 1403.0702], 

[19] L. A. Pando Zayas and N. Quiroz, Left-Right Entanglement Entropy of Boundary States, 
JHEP 1501 (2015) 110, [arXiv: 1407.7057], 

[20] D. Das and S. Datta, Universal features of left-right entanglement entropy, 

arXiv:1504.0247. 

[21] T. Quella, I. Runkel, and G. M. Watts, Reflection and transmission for conformal defects, 
JHEP 0704 (2007) 095, [hep-th/0611296]. 

[22] G. Delfino, G. Mussardo, and P. Simonetti, Scattering theory and correlation functions in 
statistical models with a line of defect, Nucl.Phys. B432 (1994) 518-550, [hep-th/9409076]. 

[23] J. L. Cardy, Boundary conditions, fusion rules and the Verlinde formula, Nuclear Physics B 
324 (Oct., 1989) 581-596. 

[24] C. Bachas, 1. Brunner, and D. Roggenkamp, Fusion of Critical Defect Lines in the 2D Ising 
Model, J.Stat.Mech. 2013 (2013) P08008, [arXiv: 1303.3616]. 

[25] C. Bachas, I. Brunner, and D. Roggenkamp, A worldsheet extension of 0(d,d:Z), JHEP 1210 
(2012) 039, [arXiv: 1205.4647]. 

[26] P. Calabrese and J. L. Cardy, Entanglement entropy and quantum field theory, J.Stat.Mech. 
0406 (2004) P06002, [hep-th/0405152]. 

[27] C. Bachas and I. Brunner, Fusion of conformal interfaces, JHEP 0802 (2008) 085, 

[arXiv: 0712.0076]. 


- 25 - 


[28] V. Petkova and J. Zuber, Generalized twisted partition functions, Phys.Lett. B504 (2001) 
157-164, [hep-th/0011021]. 

[29] 1. Brunner and D. Roggenkamp, Defects and bulk perturbations of boundary Landau-Ginzburg 
orbifolds, JHEP 0804 (2008) 001, [arXiv:0712.0188], 

[30] D. Gaiotto, Domain Walls for Two-Dimensional Renormalization Group Flows, JHEP 1212 
(2012) 103, [arXiv: 1201.0767], 

[31] A. Konechny and C. Schmidt-Colinet, Entropy of conformal perturbation defects, J.Phys. 
A47 (2014), no. 48 485401, [arXiv: 1407.6444]. 

[32] G. Poghosyan and H. Poghosyan, RG domain wall for the N=1 minimal superconformal 
models, arXiv: 1412.6710. 


- 26 - 


